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Originating from the Fertile Crescent in the Middle East, barley has now been cultivated widely 
on different soil types including acid soils, where aluminium toxicity is a major limiting factor. 
Here we show that the adaptation of barley to acid soils is achieved by the modification of 
a single gene (HvAACTI) encoding a citrate transporter. We find that the primary function 
of this protein is to release citrate from the root pericycle cells to the xylem to facilitate the 
translocation of iron from roots to shoots. However, a 1-kb insertion in the upstream of the 
HvAACTI coding region occurring only in the Al-tolerant accessions, enhances its expression 
and alters the location of expression to the root tips. The altered HvAACTI has an important 
role in detoxifying aluminium by secreting citrate to the rhizosphere. Thus, the insertion of a 
1-kb sequence in the HvAACTI upstream enables barley to adapt to acidic soils. 



1 Institute of Plant Science and Resources, Okayama University, Kurashiki 710-0046, Japan. 2 Graduate School of Bioagricultural Sciences, Nagoya 
University, 464-8601, Japan. Correspondence and requests for materials should be addressed to J.F.M. (email: maj@rib.okayama-u.ac.jp). 



NATURE COMMUNICATIONS | 3:713 | DOI: 10.1038/ncomms1726 | www.nature.com/naturecommunications 

© 2012 Macmillan Publishers Limited. All rights reserved. 



1 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1Q38/ncomms1726 



Aluminium (Al) and iron (Fe) are the two most abundant met- 
als in the Earths crust. The two metals have contrasting effects 
on plant growth; aluminium is a toxic metal, whereas iron 
is essential for plant growth. Therefore, plants rooting in soils must 
employ strategies to acquire iron, but to detoxify aluminium. Plants 
have developed two distinct strategies for the acquisition of iron 
from the soil 1 . The first strategy is employed by dicotyledonous and 
non-gramineous monocotyledonous species, and is characterized 
by the release of proton (H + ) and reductants/chelators (for exam- 
ple, electron (e), organic acids and phenolics) into the rhizosphere, 
thereby increasing iron solubility 1 . By contrast, Fe- acquisition 
by the second strategy, which is employed only by gramineous spe- 
cies, involves the secretion of hexadentate Fe 3 + -chelating substances 
(phytosiderophores as mugineic acids) and subsequent uptake of 
Fe 3 + chelates 2 . Most of the genes responsible for the different steps 
of iron acquisition have been isolated 1 . 

In contrast to Fe, aluminium is not essential to plants. Although 
metallic aluminium and aluminosilicate minerals are nontoxic, 
ionic aluminium (mainly Al 3 + ), which is released into the soil 
under acidic conditions, is toxic to all living cells. This form of Al 
rapidly inhibits root elongation at the micromolar level by targeting 
multiple cellular sites, including the cell wall and the plasma mem- 
brane, and various cellular processes, such as signal transduction 
pathways and calcium homeostasis 3 " 4 . Consequently, aluminium 
inhibits the uptake of water and nutrients and increases the sensitiv- 
ity of plants to various stresses, especially drought. Aluminium tox- 
icity has, therefore, been recognized as a major factor limiting crop 
production on acid soils, which comprise 30-40% of arable soils in 
the world 5 . However, some plant species or accessions have evolved 
strategies to cope with Al toxicity 6-8 . One of the most common 
mechanisms is the secretion of organic acid anions, such as citrate, 
oxalate and malate, from the roots 9-11 . These organic acid anions 
are able to chelate Al in the rhizosphere, thus preventing Al from 
binding to the root cell components 9 . Many plant species including 
both dicots and monocots are known to adopt this mechanism 9-11 . 
Genes encoding transporters for the Al-induced secretion of malate 
and citrate have been identified. The malate transporter gene 
ALMT1 was isolated in wheat and Arabidopsis 12 ' 13 . Genes respon- 
sible for the Al-induced citrate secretion have also been identified in 
barley, sorghum, and other plant species 14-19 . 

Barley (Hordeum vulgare), the fourth most important cereal crop 
for human food and animal feed, is one of the most Al-sensitive spe- 
cies among the small grain cereals; however, there is a wide variation 
in Al tolerance among barley cultivars 20 . This variation has been 
found to correlate with the amount of citrate secretion from the 
roots 20 . Furthermore, the Al-induced secretion of citrate from bar- 
ley roots is mediated by HvAACTl (Al-activated citrate transporter 
1) 14 , which is a member of the multidrug and toxic compound extru- 
sion (MATE) family. HvAACTl is localized to the plasma membrane 
of the root epidermal cells 14 . Although the expression of HvAACTl 
is not induced by Al, its expression level is constitutively higher in 
the Al-tolerant cultivars consistent with greater citrate secretion 14 . 
However, the mechanisms underlying the different expressions of 
HvAACTl in different barley cultivars are unknown. 

Here we report that a 1-kb insertion in the upstream enhances the 
expression of HvAACTl in the Al-tolerant cultivars. Furthermore, we 
found that this insertion also alters the spatial expression of HvAACTl 
from the original location in the mature root zone (required for 
the root-to-shoot translocation of iron) to the root tips, which are the 
sensitive targets of Al toxicity. Evolutional analysis showed that the 
insertion occurs only in the Al-tolerant cultivars from Japan, Korea 
and China, a result of adaptation to the acid soils in these areas. 

Results 

Comparison of the HvAACTl genome sequence. We have 
previously shown that the Al-sensitive barley cultivar Morex and 



the Al-tolerant cultivar Murasakimochi differed in the level of 
HcAACTl expression 14 . To investigate the underlying mechanism, 
we analysed the genome sequences of HvAACTl in both cultivars. 
We found an insertion with 1,023 bp in the upstream of HvAACTl 
coding region in Murasakimochi, but not in Morex (Fig. la, 
Supplementary Fig. SI). HvAACTl is characterized by a large 
intron (about 6kb) in the 5'-untranslated region (UTR) region and 
this insertion is localized at 4.8 kb upstream from initiating codon 
(ATG) in the Murasakimochi genome (Fig. la, Supplementary Fig. 
SI). We then checked this insertion in other cultivars and found 
that it is present in all Al-tolerant cultivars tested (Fig. lb). Database 
search revealed that the insertion shows similarity with the CACTA- 
like transposon in rice. Furthermore, this insertion has a terminal 
inverted repeats (17 bp) at the 5' end and 3' end (Fig. lc). 

Determination of the transcriptional start sites. We determined 
the potential transcriptional start sites of HvAACTl using 5'-RACE. 
Although there were three different 5'- UTR sequences in Morex, 
there was only one transcriptional start site (TSS) at 6kb upstream 
from ATG of HvAACTl (Fig. ld,e). By contrast, there were multi- 
ple TSSs in Murasakimochi in addition to the same TSS found in 
Morex. These TSSs originated from the insertion (Fig. lc), yield- 
ing various mRNAs with different 5'-UTR during splicing processes 
(Fig. ld,e). In addition, these 5'-UTR sequences did not contain an 
initiating codon. All transcripts encode the same functional gene 
(Fig. le). 

Promoter activity assay. To investigate the function of this inser- 
tion, we determined the promoter activity by transiently introduc- 
ing green fluorescent protein (GFP) fused with a part (262 bp) of the 
insertion before most TSSs into onion epidermal cells. A putative 
common promoter and the 35S promoter were also fused with GFP 
for comparison. GFP signal was observed in the cells introducing 
all promoters (Fig. If). These results indicate that the insertion 
functions as a promoter in onion epidermal cells. 

Expression patterns. To investigate how this additional promoter 
in the insertion region affects the expression of HvAACTl, we per- 
formed a detailed expression analysis of the 5'-UTR region. The 
HvAACTl expression level gradually increased behind the multi- 
ple TSS sites (Fig. 2). Spatial analysis showed that the expression 
level at the common UTR was very low in the root tips (0-3 mm) 
of both Murasakimochi and Morex (Fig. 3a-c). In the other parts 
of the roots, the expression level at the common UTR in Morex was 
slightly higher than or similar to that of Murasakimochi. At the 
insertion-specific region, the expression was detected only in Mura- 
sakimochi, with both the root tips (0-3 mm) and the other parts of 
roots showing similar levels of expression (Fig. 3c). For the coding 
region, the expression level was 50 times higher in Murasakimochi 
than in Morex in the root tips and 3 times higher in the other parts 
of roots (Fig. 3c). 

Tissue-specific expression of HvAACTl. To examine the tissue- 
specific expression of HvAACTl, we obtained RNA from the central 
cylinder (pericycle and inner tissues) and the outer tissues (cor- 
tex and epidermis) of root tissues by using laser microdissection 
(Fig. 3d-f ). In the root tips (around 2 mm from the apex), the expres- 
sion of HvAACTl was hardly detected in either the central cylin- 
der or the outer tissues of Morex (Fig. 3g). By contrast, HvAACTl 
was expressed at a higher level in both tissues in Murasakimochi 
(Fig. 3g). In the mature zone (around 20 mm from the apex) where 
xylem is well developed, the expression level of HvAACTl in the 
central cylinder was similar between Morex and Murasakimochi 
(Fig. 3g). However, in the outer tissues, no expression of HvAACTl 
was detected in Morex, whereas a high expression was found in 
Murasakimochi. 
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Figure 1 1 Location and function of the upstream insertion, (a) Schematic presentation of the genomic structure of HvAACTI in an Al-sensitive cultivar 
(Morex; Mo) and an Al-tolerant cultivar (Murasakimochi; Mu). White and colour vertical boxes show UTR exons. Black vertical boxes and black lines 
show ORF-exons and introns, respectively. Grey line shows a 1,023-bp insertion at the 4.8 kb upstream of the coding region in Murasakimochi. TSS, 
transcriptional start site. Grey arrows indicate the positions of PCR primers used for the insertion detection in b. (b) Detection of the 1-kb upstream 
insertion in different barley cultivars. 1, Morex; 2, Golden Promise; 3, BC95; 4, Z504; 5, K0046; 6, BC282; 7, BC366; 8, ALP3; 9, BC26; 10, Murasakimochi; 
11, BC262; 12, BC191; 13, Akashinriki. (c) Sequence of HvAACTI upstream insertion. White triangles show multiple HvAACTI TSS and the number means 
frequency of identification of TSS by cloning. Sequence of blue character indicates 5'-UTR-2 in Murasakimochi. The insertion region contains a terminal 
inverted repeat (TIR) of 17 bp at the 5' and 3' ends, (d) Image of agarose gel showing the 5'-RACE PCR products. HvAACTI gene-specific primer and 
RNA oligo primer were used for amplifying HvAACTI 5'-UTR region, (e) Schematic presentation of 5'-UTR variants of HvAACTI between Morex and 
Murasakimochi. The sequences were read from 20 5'-RACE clones. White and colour boxes show the 5'-UTR and colours corresponding to diagram in 
a. The number in the box shows 5'-UTR position on the genome. Box 2' includes variations from different transcriptional start sites, (f) Promoter assay. 
Common promoter (PI in a), 1-kb insertion-originated region (P2 in a), or 35S promoter (as a control) was fused to GFP and then transiently introduced 
to onion epidermal cells. Scale bars, 100 urn. 



Transformation of HvAACTI into barley. To further investigate the 
role of the 1-kb insertion in Al tolerance, we introduced HvAACTI 
into an Al-sensitive cultivar (Golden Promise) under the control of 
three different promoters; the specific insertion region (B3), pro- 
moters of Murasakimochi (Bl) or Morex (B2) (Supplementary 
Fig. S2). Analysis with TO plants showed that there was a positive 
correlation between the HvAACTI expression level and citrate secre- 
tion in transgenic plants carrying the Murasakimochi promoter 
(Bl) and the specific insertion region (B3; Supplementary Fig. S3). 
We randomly selected four independent lines for further analysis. 
In the Tl transgenic lines, the expression level of HvAACTI in the 
root tips (0-3 mm) was significantly enhanced in the lines carrying 
HvAACTI driven by the Murasakimochi promoter or the specific 
insertion region (Fig. 4a,d), but was unaffected in those carrying 
the Morex promoter. Similarly, Al-induced secretion of citrate and 
Al tolerance were significantly enhanced in the lines driven by the 
Murasakimochi promoter or the insertion region (Fig. 4b,c,e,f), but 
not in those driven by the Morex promoter. These results indicate 
that the specific insertion enhances the expression of HvAACTI in 
the root tips, resulting in enhanced Al-induced secretion of citrate 
and subsequently Al tolerance. 



Role of HvAACTI in iron translocation. The similar expression level 
of HvAACTI in the central cylinder tissues of the mature root zone of 
both cultivars suggests that HvAACTI has other functions (Fig. 3g). 
HvAACTI is a homologue of OsFRDLl in rice and AtFRD3 in Ara- 
bidopsis 21 ' 22 , which are localized in the pericycle cells and involved 
in the translocation of iron from the roots to the shoots by transport- 
ing citrate into the xylem. To test whether HvAACTI is also involved 
in iron translocation, we introduced the Morex-type HvAACTI into 
the rice mutant osfrdll defective in the translocation of iron 21 . When 
the mutant was grown at a low Fe concentration (0.2 |iM), the new 
leaves of osfrdll mutant showed chlorosis typical of Fe deficiency 
(Fig. 5a), whereas the transgenic lines carrying HvAACTI rescued 
the phenotype of the mutant. The SPAD (indicator of chlorophyll 
content) value of the youngest expanded leaf was lower in the osfrdll 
mutant than the wild-type rice; this was restored to the wild-type 
level by the transformation with HvAACTI (Fig. 5b). In the presence 
of high Fe (10jaM), all lines showed similar SPAD values (Fig. 5c). 

In the transgenic lines, the concentrations of citrate and iron in 
the xylem sap were significantly higher than those of the osfrdll 
mutant and were comparable to those of the wild-type rice (Fig. 5d,e). 
Expression of HvAACTI also increased the iron concentration 
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Figure 2 | HvAACTI expression level at different 5'-UTR regions. 

(a) A schematic model of primer position and gene structure. The colour 
arrows show regions for PCR amplification in b. Black arrows show TSS. 

(b) HvAACTI expression at different 5'-UTR region. The colour of column 
corresponds to that in a. Statistical comparison was performed by one- 
way ANOVA followed by Dunnett's multiple comparison test. All data 
were compared with the expression of Morex ORF in the -Al treatment. 
*P< 0.05. Data are shown as the mean values + s.d. (n = 3). 



in the shoots, but decreased the iron concentration in the roots 
(Fig. 5f,g). There were Fe precipitates in the osfrdll mutant roots; 
these disappeared in the transgenic lines (Fig. 5h). These results 
indicate that HvAACTI also has a role in transporting Fe from the 
roots to the shoots by releasing citrate to the xylem. 

Evolutional analysis of HvAACTI. To examine the evolutionary 
origin of the 1-kb upstream insertion, we investigated its presence 
in 154 and 265 accessions of wild barley (Hordeum vulgare ssp. 
spontaneum) and cultivated barley {Hordeum vulgare ssp. vulgare), 
respectively, which have been collected from different geographic 
locations around the world. The insertion was not found in any of the 
wild barley accessions examined (Table 1). It was detected only in 20 
cultivated barley accessions distributed in the acid soil areas in Japan, 
Korea and China (Table 1, Supplementary Fig. S4). We compared 
the Al tolerance of these accessions with 20 other accessions ran- 
domly selected from the collection by Eriochrome cyanine staining 
method, which reflects Al accumulation in the roots (Supplementary 
Fig. S5). All cultivars with the insertion showed higher Al tolerance 
(less Al staining) compared with those without the insertion. 

In order to address the origin and the spread of the 1-kb inser- 
tion in the global barley population in more detail, we performed 
a haploytpe analysis by sequencing the region spanning 1,237 bp 
within the HvAACTI gene in 48 barley accessions (Supplementary 
Table SI). There were only nine single -nucleotide polymorphism 
and three Indel sites (including the 1-kb insertion) in this region 
among the accessions tested (Supplementary Table S2). Ten haplo- 
types can be defined based on the polymorphisms (Supplementary 
Table S2). Based on the silent base substitutions, we estimated that 
the nucleotide diversity, 71, of the silent sites in the HvAACTI gene 
to be 0.0024 and 0.0007 in wild and cultivated barley, respectively. 
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Figure 3 | Altered expression patterns of HvAACTI by a 1-kb upstream insertion, (a) Schematic presentation of sampling position of roots, (b) Primer location 
designed for spatial expression analysis. The white and colour boxes showed 5'-UTR sequences of Morex (Mo) and Murasakimochi (Mu). Small black arrows 
show TSS. Large black arrows show regions for PCR amplification. The number in the boxes shows the position of 5'-UTR. (c) Spatial expression of HvAACTI in 
different root zones assayed with different primer sets at the common 5'-UTR, the insertion region and the coding region of HvAACTI, respectively, as shown 
in b. (d) Schematic presentation of different tissue samples separated by laser microdissection. The epidermis and cortex were collected as outer tissues 
and inner tissues from pericycle as central cylinder. (e,f) Outer tissues (e) and central cylinder (f) of the roots after laser microdissection, (g) Tissue-specific 
expression of HvAACTI. Expression level relative to 2 mm central cylinder of Morex is shown. Data in c and g are shown as the mean values+s.d. (n = 3). 
Statistical comparison was performed by one-way ANOVA followed by Dunnett's multiple comparison test, (c) All data were compared with 0-3 mm Morex 
assayed with the ORF primer. *P< 0.05. (g) The expression was compared with the central cylinder of Morex at 20 mm. *P< 0.05. 
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Figure 4 | Function of the upstream insertion in transgenic barley, (a) Expression of HvAACTl in the root tips (0-3 mm). Plasmids containing HvAACTl 
ORF fused with Murasakimochi promoter (B1), Morex promoter (B2) or insertion-originated promoter (B3) were transformed into the barley cultivar 
Golden Promise (GP). (b) Al-induced secretion of citrate, (c) Al tolerance in transgenic barley lines, (d-f) Expression of HvAACTl in the root tips 
(0-3 mm) (d), Al-induced secretion of citrate (e) and Al tolerance in transgenic barley lines (f) in a separate experiment using different independent 
lines. Data are shown as the means + s.d. (n = 3) for a bde and (n = 10) for c,f. Statistical comparison was performed by one-way ANOVA followed by 
the Tukey's multiple comparison test. All data were compared with GP. *P< 0.05. 



A haplotype network was constructed to examine whether the 
1-kb insertion is of single or multiple origins (Fig. 6). The acces- 
sions with the 1-kb insertion belong to the haplotype 10, which 
originated from the haplotype 1 (Supplementary Fig. S6). The rela- 
tionship between haplotype 10 (with insertion) and haplotype 1 
(without insertion) is also supported by the phylogeny of substitu- 
tion mutations (Supplementary Fig. S6) derived from the Supple- 
mentary Table S2. These results indicate that the 1-kb insertion is 
from a single origin during the expansion of barley cultivation onto 
the acid soil areas. 

Discussion 

Barley was domesticated from its wild ancestral form H. vulgare 
ssp. sponatneum in the Near-East about 10,000 years ago 23 . Soils 
in this region are predominantly alkaline, where there is no threat 
of Al toxicity stress, whereas the Fe availability is limited. Therefore, 
barley has developed an efficient strategy for acquiring insoluble 
iron in alkaline soil by synthesizing and secreting large amount of 
phytosiderophores 2 . Around 3,000 years ago barley cultivation was 
expanded to East Asia, where the crop started to encounter Al tox- 
icity stress on acid soils. Our results indicate that cultivated barley 
acquired Al tolerance through a 1-kb insertion in the upstream of 
an Al-tolerant gene, HvAACTl. 

HvAACTl encodes a citrate transporter 14 . Its original role in 
barley is to release citrate to the xylem for efficient translocation 
of Fe (Fig. 7). Citrate forms complexes with iron in the xylem to 
prevent Fe precipitation 21-22 . The release of citrate to the xylem is 
mediated by citrate transporters such as AtFRD3 in Arabidopsis 22 
and OsFRDLl in rice 21 ; the barley HvAACTl is also a homologue 
of these citrate transporters. They are expressed in the pericycle 
cells in the mature root zones, where xylem is well developed. 



When HvAACTl was introduced into the rice mutant osfrdll, it 
restored all mutant phenotypes, including the concentrations of 
citrate and Fe in the xylem sap, the leaf chlorophyll content and Fe 
precipitation in the roots, to the wild-type levels (Fig. 5). Therefore, 
HvAACTl in the root pericycle cells functions as a citrate trans- 
porter required for the Fe nutrition; this role is conserved in all 
barley accessions regardless of their Al tolerance (Figs 3 and 7). In 
fact, HvAACTl also showed a transport activity for citrate in the 
absence of Al although the activity was greatly enhanced by Al 
exposure 14 . 

Interestingly, HvAACTl gains a second function with the 1-kb 
insertion in the upstream of the coding region in the Al-tolerant 
barley cultivars (Fig. 7). This insertion functions as a promoter 
because of the presence of multiple transcriptional start sites 
(Figs 1 and 2). The insertion not only enhances the expression of 
HvAACTl in the whole root, but also spatially alters its expression 
pattern by markedly increasing the expression in the root tips of 
the Al-tolerant cultivars (Fig. 3). Root tips are the most sensitive 
target of Al toxicity 24 . The enhanced expression of HvAACTl in this 
root zone results in increased citrate secretion to the rhizosphere, 
thereby protecting the root tips from Al toxicity (Fig. 4). 

Gene expression is regulated by different mechanisms, such as 
chromatin condensation, DNA methylation, transcriptional ini- 
tiation, alternative splicing of RNA, mRNA stability and so on 25 . 
Among the Al-tolerant genes identified, a minimum MITE insertion 
size in the promoter region of SbMATE in sorghum was suggested to 
be involved in regulating the expression of this gene 15 . On the other 
hand, tandemly repeated elements of length 33-803 bp and located 
in the upstream of the TaALMTl -coding region are associated with 
a high expression of TaALMTl in wheat 26-27 . The mechanism for 
regulating HvAACTl expression in barley appears to be distinct 
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Figure 5 | Complementation test for iron translocation in rice, (a) Rescue of leaf chlorosis by the introduction of Morex HvAACTl (b,c) SPAD value of 
the youngest expanded leaf of transgenic rice carrying Morex-type HvAACTl promoter grown in a nutrient solution containing 0.2 (b) or 10 (c) jiM FeSO^ 
Data are shown as means + s.d. (n = 4-6). (d,e) Concentration of citrate (d) and iron (e) in the xylem sap of different rice lines. (f,g) Concentration of iron 
in the shoots (f) and roots (g). Data are shown as means + s.d. (n = 3-5). Statistical comparison was performed by one-way ANOVA followed by Dunnett's 
multiple comparison test. *P< 0.05. All data were compared with osfrdll. (h) Rescue of iron precipitation in the roots by the introduction of Morex 
HvAACTl Scale bars in a and h indicate 10 mm and 100 Jim, respectively. 



Table 1 1 Detection of the 1-kb insertion in HvAACTl 
upstream region in wild and cultivated barley. 

Taxon Geographical group 1 kb Insertion 



Hordeum vulgare ssp. spontaneum 



Hordeum vulgare ssp. vulgare 



Near East 


76 




Transcaucasus 


16 




Southwest Asia 


57 




Unknown 


5 




Total 


154 


0 


New continent 


5 




Europe 


42 




South Europe 


17 




North Africa 


11 




Ethiopia 


35 




Turkey 


18 




Near East 


25 




Transcaucasus 


9 




Southwest Asia 


10 




Himalaya 


23 




China 


21 


1 


Korea 


20 


5 


Japan 


9 


14 


Total 


245 


20 



from those in other plant species. Furthermore, TaALMTl and 
SbMATE are responsible for Al tolerance only, whereas HvAACTl 
has evolved dual functions; Fe translocation and Al tolerance. 

6 



The cultivated accessions of barley are genetically divided into 
the Western and Eastern groups 28 . The occurrence of the 1-kb inser- 
tion in the upstream region of HvAACTl gene might originate in 
the Eastern region where acid soils are distributed (Table 1). This is 
supported by the low nucleotide diversity, 71, of HvAACTl, which 
is markedly lower than those reported previously for other genes 28 . 
The low nucleotide diversity of HvAACTl gene indicates that it is 
functionally and evolutionally conserved in barley species and that 
the 1-kb insertion happened during the expansion of barley cultiva- 
tion into the acid soil areas. This is reasonable considering that the 
original role of HvAACTl is to facilitate the translocation of iron. Our 
study provides evidence for the acquisition of Al tolerance as one of 
the adaptive traits; it is this adaptation that allows the expansion of 
barley cultivation onto the acidic soils in the Far-Eastern region. 

Methods 

Plant materials and growth conditions. Barley (Hordeum vulgare) seeds were 
soaked in deionized water overnight at 20 °C in the dark and then transferred 
to nets floating on a 0.5 mM CaCl 2 solution in a 1.5-liter plastic container. After 
growth at 20 °C for 3-4 days, the seedlings were transferred to a 3.5-liter plastic 
pot containing one-fifth Hoagland solution. 

Genomic sequence analysis. A BAC clone containing HvAACTl was screened 
by using primers: 5'-TGGAGGAAGCATAGTATC-3' and 5'-CACCTGGAGGT 
ATGAA-3'. The shotgun library of the selected BAC was sequenced and assembled 
to obtain the genomic DNA sequence of HvAACTl . 

The upstream of the HvAACTl genome sequence was obtained from an 
Al-sensitive cultivar (Morex) and an Al-tolerant cultivar (Murasakimochi) by 
PCR using primers; 5 , -AACAGCTCTGCACCCAGTCGCCGGCCACTT-3 , 
and 5'- GATGTTC ATG AC AAAGAGGTAG AG- 3'. The fragments obtained were 
directly sequenced using BigDye Terminators V3.1 cycle sequencing kit on an 
ABI PRISM 310 genetic analyser. 

To examine the presence of the insertion in other accessions, 13 accessions dif- 
fering in Al tolerance, including Morex and Murasakimochi, were used. Fragments 
including the insertion were amplified by PCR using Takara Ex Taq (Takara) from 
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Figure 6 | Haplotype network among ten haplotypes of HvAACTI. 

Haplotype network is constructed using both substitution and Indel 
polymorphisms by the maximum-parsimony criterion (heuristic search) in 
PAUP*. Thin and bold lines represent base substitution and Indel mutations, 
respectively. Asterisks indicate inferred intermediate haplotypes. Numbers 
above or next to the lines indicate the variable position. A 1-kb insertion 
in the 5'-flanking region is marked by an arrow. One homoplasious 1-bp 
insertion is marked by dotted lines. Observed haplotypes in wild barley are 
shown in circles, and in cultivated barley in squares. 



the genomic DNA extracted from each cultivar. The primers used were 5'-GGTC 
CAACACTCTACCCTTCCTT-3 , and 5'-GGTGCGAGTTGCCCCTAGCTATTA 
CAGA-3'. The PCR products were separated by electrophoresis on a 1.0% agarose 
gel using 0.5x TBE buffer stained with ethidium bromide. 

Identification of the TSS. To determine the TSS of HvAACTI, total RNA from 
the roots of both Morex and Murasakimochi was isolated with an RNeasy plant 
Mini Kit (Qiagen). 5'-RACE was performed by GeneRacer kit (Invitrogen) using 
2 jig total RNA. The ligated mRNA with RNA oligo adapter was transferred to a 
first-strand cDNA using Superscript III RT and the random primers in the reverse- 
transcription reaction. The 5' ends of HvAACTI cDNA were amplified by PCR 
with the GeneRacer 5' primer and gene-specific primer 5'-CTCTACCTCTTT 
GTCATGAACATC-3', GeneRacer 5' Nested primer and gene-specific primer 
5'-GATGTTCATGACAAAGAGGTAGAG-3'. The PCR products were cloned into 
pCR2.1 TOPO vector (Invitrogen) and cloned plasmids were sequenced. 

Transient promoter activity assay. To investigate the promoter activity two dif- 
ferent regions were fused with GFP. Putative promoter 1 was located at 1 kb from 
just after the common TSS, which was amplified by KOD-Plus (Toyobo) using 
the forward primer, 5'-ACGGACACCTTCGGAGATACCTGTAGAGCA-3' and 
reverse primer, 5'-CGCCTATTTATAGGGGCCACACGACCCACCATGG-3'. 
Putative promoter 2 (262 bp) was located in the specific insertion region just after 
the most TSS and was amplified using the forward primer, 5'-ATATGCCGT 
GGACAGA ATAC AGTTATGCCC- 3' and reverse primer, 5'-CCATGGCCT 
GGAGTAGGTTGTCTCTCG ATGC AGCG- 3'. The fragment amplified was 
introduced into pGEM vector by pGEM-T (Easy) Vector Systems (Promega). These 
fragments were introduced into the GFP-nopaline synthase (NOS) terminator plas- 
mid, which has the native promoter deleted. As a positive control, GFP was fused 
to the 35S promoter. The constructs were then bombarded into onion epidermal 
cells. After overnight incubation, GFP fluorescence was observed with a confocal 
laser scanning microscopy (LSM700; Carl Zeiss). 

Expression pattern analysis. Roots of both Morex and Murasakimochi grown in 
a 0.5-mM CaCl 2 solution for 4 days were separated into different segments (0-3, 
3-6, 6-9, 9-12, 12-20 mm) with a razor. The samples were ground in liquid nitro- 
gen and total RNA was isolated with RNeasy plant Mini Kit (Qiagen). cDNA was 
synthesized as described above. The expression levels of HvAACTI in different root 
zones were quantified by absolute quantitative real-time PCR using different prim- 
ers. For the common UTR region, primers used were 5'-GCCAGTAGCTCTCT 
GCGCTC-3' and 5'-CCTGAAGTGGCGAAGTGG-3'. For the specific insertion 
region, primers used were 5'-GCTGCATCGAGAGACAACCT-3' and 5'- TCAGA 
AGCTGTCTGCTCGGG-3', and for the coding region, the following primers were 
used: 5'-TTCGTCGCTGAAGAAGATGC-3' and 5'-GCACATTGCAGGCATCT 
GAA-3'. Reactions were run on the MasterCycler ep realplex (Eppendorf). To 
generate standard curves for absolute quantification for each region, a series of 
dilutions (from 1x10" 1 to 10 _6 ng) of plasmids were made and then assayed by 
real-time PCR. Cj values for each sample were converted into absolute copy 
numbers using the standard curves. 



Morex genome / 
(Al sensitive) / ' 



Murasakimochi genome 
(Al tolerance) 




ii Minimi 

ii 



Jjo-3 mm 

Fe translocation 



J 0-3 mm 

Al tolerance 



Figure 7 | Schematic presentation of HvAACTI roles in barley. HvAACTI 
in the pericycle cells of root mature zone releases citrate to the xylem, 
which is required for the translocation of Fe from the roots to the shoots in 
both Al-sensitive and -tolerant cultivars. A 1-kb insertion in the upstream in 
the Al-tolerant cultivars alters the expression to the root tips and enhances 
the expression in both the root tips and mature zones. HvAACTI expressed 
in the root tips releases citrate to the rhizosphere to detoxify Al. 

To investigate tissue-specific expression of HvAACTI, laser microdissection 
was performed by the method of Takahashi et al. 29 Roots of both Morex and 
Murasakimochi were dissected at 0.5-3.5 mm and 18.5-21.5 mm to collect antero- 
posterior 1.5 mm in a surplus. The samples were immediately immersed in a fixing 
solution containing ethanol/acetic acid at a 3:2 ratio. Central cylinder and outer 
tissues were collected from the root tissue sections using a Veritas Laser Microdis- 
section System LCC1704 (Molecular Devices) and used for total RNA extraction. 
The cDNA from each sample was used for the quantitative analysis of gene 
expression. Cyclophilin was used as an internal control by the AACt method. 
HvAACTI expression levels were normalized with that in the 0-3 mm central 
cylinder of Morex. 

Preparation of transgenic barley and rice. The promoter region was amplified 
by 5'-GGTACCACGGACACCTTCGGAGATACCT-3' and 5'-GGACGTCGAC 
GACGGCCACCCACTTCTTGT-3' primers for the Bl fragment, 5'-TACCT 
GCAGGCAGTAGCTCTCTGCGCTCAG-3' and 5'-GGACGTCGACGACG 
GCCACCCACTTCTTGT-3' primers for the B3 fragment from Murasakimochi 
genomic DNA, using KOD-FX (Toyobo; Supplementary Figs SI and S6). The B2 
promoter region was amplified by 5'-GGTACCACGGACACCTTCGGAGA 
TACCT-3' and 5'-ACAAGAAGTGGGTGGCCGTCGTCGACGTCC-3' primers 
from the Morex genomic DNA (Supplementary Figs SI and S6). These fragments 
were cloned into pCR2 vector by TOPO XL PCR Kit (Invitrogen). The HvAACTI 
open reading frame (ORF) was amplified by using forward primer 5'-ACTCCTCT 
GTCCGAGCACGA-3' and reverse primer 5'-CCCGGGCTCCTCGCACAAAAG 
GCAGACGACC-3' by using KOD-FX from the Morex cDNA. This fragment 
was cloned into pGEM vector (Promega) and sequenced. The ORF was inserted 
between with NOS terminator into pUC19 vector. The plasmids of cloned 
promoters as described above were digested by Kpnl and BamHl and the plasmid 
of cloned ORF-NOS terminator was digested by Notl and BamHl. These fragments 
were ligated into pPZP vector 30 and then transformed into Agrobacterium 
tumefaciens (strain EHA101). Immature embryos of barley cultivar Golden 
Promise were used for Agrobacterium-mediated transformation 31-33 . We 
obtained several independent transgenic lines and three or four each were 
randomly selected for further analysis. Transformation of the Morex promoters 
fused with HvAACTI cDNA into a rice mutant osfrdll defective in iron 
translocation 21 was performed 34 . 

Phenotypic analysis of transgenic barley and rice. To measure Al tolerance, 
transgenic barley lines were exposed to a 1-mM CaCl 2 solution with or without 
3 \lM Al at pH 5.0. Root elongation during 24 h was measured with a ruler. Root 
exudates of the transgenic lines and WT (3 weeks old) were collected by exposing 
plants to a 1-mM CaCl 2 solution containing lOuJVl Al for 6h. Citrate was meas- 
ured by an enzymatic method 35 . 

Transgenic rice lines, WT and osfrdll mutant were grown in a one-half strength 
Kimura B solution containing 0.2 uM or 10 uM FeS0 4 at 25 °C. The SPAD value 
of the youngest fully expanded leaf was measured after growth for 15 days. Plants 
grown in the presence of 0.2 |iM FeS0 4 were decapitated at 2 cm above the roots 
and the xylem sap was collected from the cut end for 1 h by micropipette. The 
concentrations of citrate and iron in the xylem sap were determined as described 
previously 21 . The Fe concentrations in the roots and shoots were determined by 
ICP-MS (Aglient Technologies). Perls' staining was performed in plants grown in 
a solution containing 10 [iM FeS04 for 1 month 21 . 
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Geographic distribution analysis. A total of 265 cultivated barley accessions 
collected from around the world and 154 wild barley accessions were used to 
examine the presence of the specific insertion. Primers for amplified fragments 
including the insertion region were used. Genomic DNA was extracted from each 
accession and used as a template. Touchdown PCR was run by using KOD-FX, 
followed by standard PCR performed after 30 cycles. 

Al tolerance was evaluated by Eriochrome cyanine staining 36 . A total of 40 
cultivars including 20 with the insertion and 20 without the insertion, which were 
randomly selected from the collection, were used. Seedlings (1 -week-old) were 
exposed to a 1.0-mM CaCl 2 solution (pH 5.0) for 48 h and then stained in 0.1% 
Eriochrome cyanine for 15 min 36 . Roots were then washed with deionized water 
and photographed. 

The distribution of barley accessions with or without the insertion was mapped 
onto a world soil atlas with different soil pHs (http://www.sage.wisc.edu/atlas/ 
maps.php?datasetid=20&includerelatedlinks= l&dataset=20). The information 
of the accession origin of cultivated barley was based on Barley DB (http://www. 
shigen.nig.ac.jp/barley/collectionMap/map.html). Wild barley accessions have 
been described by Konishi 37 (Supplementary Table S3). 

Phylogenic analysis. To generate a neighbour-joining tree, we sequenced the 
region spanning 1,237 bp within the HvAACTl gene (688 bp from flanking region 
to the insertion and 549 bp from the second intron) in 48 barley accessions 
comprising of 40 cultivated ones and 8 wild ancestors, H. spontaneum (Supplemen- 
tary Table SI). Alignment of all sequences was performed by ClustalW. Regions 
containing gaps and missing data were not included. Neighbour-joining tree was 
constructed with MEGA5. The evolutionary distances were computed using the 
Kimura 2-parameter method. Bootstrap test with 1,000 replicates was also carried 
out by MEGA5. 

Haplotype network analysis. Haplotype network was constructed using both 
substitution and Indel polymorphisms by the maximum-parsimony criterion 
(heuristic search) in PAUP* (ref. 38). Indel mutations were recoded and treated 
as a fifth base, and all the Indel sites were considered as a single mutation. 
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